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Abstract 

The transition from the model of a long Josephson junction of variable width to 
the model of a junction with a coordinate-dependent Josephson current amplitude is 
effected through a coordinate transformation. This establishes the correspondence 
between the classes of Josephson junctions of variable width and quasi-one-dimensio- 
nal junctions with a variable thickness of the barrier layer. It is shown that for a 
junction of exponentially varying width the barrier layer of the equivalent quasi- 
one-dimensional junction has a distributed resistive inhomogeneity that acts as an 
attractor for magnetic flux vortices. The curve of the critical current versus magnetic 
field for a Josephson junction with a resistive microinhomogeneity is constructed 
with the aid of a numerical simulation, and a comparison is made with the critical 
curve of a junction of exponentially varying width. The possibility of replacing a 
distributed inhomogeneity in a Josephson junction by a local inhomogeneity at the 
end of the junction is thereby demonstrated; this can have certain advantages from 
a technological point of view. 

1. INTRODUCTION 

An interesting research topic in the physics of Josephson junctions (JJs) and one that 
has received intensive development in recent years is the influence of a potential due to 
an external field or to the geometry on the motion of vortices in long junctions. As 
examples one can cite the study of the potential created by a spatial variation of the 



*E-mail: elis@jinr.ru 
'I'E-mail: todorlb@jinr.ru 
■l-E-mail: shukrinv@thsunl.jinr.ru 



radial component of the magnetic field in an annular junction [1 and of a field-induced 
saw-tooth potential in the same annular junction (2J and the study of the fluxon qubit in 
a modified annular (heart-shaped) geometry jS] . Potentials produced by spatial variation 
of the width of the junction, corresponding to a force acting on the fluxon in the direction 
of narrowing, were studied in Refs. ^ - ^Tj. Unidirectional motion of fiuxons has a 
certain advantage in the design of oscillators based on the motion of fiuxons (121 El E]- 
In particular, an exponential variation of the junction width provides better impedance 
matching with an output load [4^ and allows one to avoid the chaotic regimes inherent to 
rectangular junctions. The structure and stability of static magnetic flux vortices and the 
corresponding critical curves in long JJs of exponentially varying width were investigated 
in Refs. 1113110], and 

In this paper it is shown that effects arising in a junction of variable width also 
take place for JJs whose barrier layer contains a resistive inhomogeneity. The critical 
curves — critical current versus magnetic field — are constructed numerically for long 
JJs of exponentially varying width and rectangular junctions with a localized resistive 
inhomogeneity. From a technological standpoint JJs with a resistive inhomogeneity can 
be preferable in some cases. 



2. COORDINATE TRANSFORMATION 

We consider a long JJ whose dimension along the y axis (width) is a smooth function 
W{x) on the interval x G [0, L]. Here L is the length of the junction (all of the quantities 
used are rendered dimensionless in the standard way [T3j). We assume that 

meLxW{x) ^ L . 

x£[0,L] 

We introduce the shape function a{x) of the junction with the aid of the relation 

"^") = -W(^' 

Here and below a subscript is used to denote differentiation with respect to the correspon- 
ding independent variable. The case Wx{x) = corresponds to a rectangular junction. 
Then the boundary- value problem for the static magnetic flux (p{x) in such a J J is written 
in the form 



-<^xx + (t{x) + simp + g{x) = , (2a) 
<Px{0) = ho, iPx{L) = hi . (2b) 

The concrete expressions for q{x), ho, and hL depend, in particular, on the means of 
injection of the external current 7. For junctions with an overlapping geometry, as are 
considered below, the base current 7(2;) flows through the whole junction. Consequently, 
g{x) = 7(x) —ahs, ho = hs, hi = hs {Jib is the external magnetic field, which is directed 
along the y axis in the plane of the barrier layer). The term —ahs is due to the variable 
width of the junction. Neglecting surface effects, one can to a first approximation set 
7(0;) = const. A derivation of Eq. (Paj) in the general case is given in Ref. j^]. 
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Let us consider an ordinary change of spatial coordinate: 



e = v{x) , (3) 

where the variable ^ G [0,/], and the right-hand side is a solution of the boundary- value 
problem 

-Vxx + (j{x)vx = , (4a) 

v{0) = 0, min = 1 . (4b) 

xe[o,L] 

The parameters L and I are connected by the relation I = v{L). The first boundary 
condition in (4b) matches up the initial points of the intervals of variation of the variables 
X and ^, while the second serves for normalization (see below). 

In the new independent variable the boundary value problem Q takes the form 

-</'€5+Jz)(0sin¥? + j(0 = 0, (5a) 
^^{0)=Ho, ^^{1)=Hl. (5b) 

Here joiO = "^1(0? i(0 = ""1(0 ""(0 is the inverse function of v{x), and also 

Ho = u^{0)ho, HL = u^{l)hL. 



From a formal standpoint Eq. ()5aj) describes the magnetic flux distribution in a 
one-dimensional J J with a variable amplitude joiO of Josephson current. Such a 
junction is inhomogeneous [121 — thickness of its barrier layer varies along the length 
as d{^) ~ M^^(^). The term j(^) is a current that varies along the junction, generated by 
variation of the geometry, external magnetic field, and external current. 

Thus the transformation 0, Q establishes a relation between a J J of variable width 
and a "one- dimensional" junction having a thickness of the insulating layer that varies 
along the junction. For brevity we shall speak of x junctions and junctions. 

For junctions of exponentially diminishing width (referred to below as EJJs), which 
were considered in Refs. jSI - [101, we have a = const > 0. Then W{x) = VFq exp(— ax), 
where Wo is the width of the JJ at the x = end. In this case the transformation (3) 
takes the form 

e = - ie'^' - 1) 



a 



Then for the right-hand boundary of the ^ junction we find / = (e — l)/a > L, i.e., the 
^ junction is longer than the corresponding x junction. 

In the particular case considered, a change of the normalized amplitude of the Joseph- 
son current in ^-junction 

j^(0 = n|(0 = — ^— 2, (6) 

(1 + ^0 

can occur as a result of a barrier-layer thickness that varies smoothly along the junction. 
Here the amplitude is maximum (and the barrier-layer thickness minimum) at the left 
end of the junction (see the second condition in Eq. ()4b|)) and minimum (maximum, 
respectively) at the right end. Thus a resistive inhomogeneity is an attractor for 
magnetic flux distributions in the junction, drawing the latter toward the right end, as 
has been observed experimentally Q and confirmed by numerical simulation |H [TUl ITT|. 
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3. MODELS OF A JUNCTION WITH A RESISTIVE 
INHOMOGENEITY 



It is of interest to investigate the possibility of replacing an inhomogeneity distributed 
over the whole length of the junction by an inhomogeneity localized near the right end. 
The simplest choice is a point inhomogeneity of "strength" > at the right end of the 
junction. In that case [111 

joix) = 1 — fi6{L — x) . 

Here S{x) is the Dirac delta function. A JJ with such an inhomogeneity will be denoted 
below as 5JJ. For /x = the junction is homogeneous. The corresponding boundary- value 
problem in the case of overlapping geometry is described in the form [x G (0, L)] 

-ifncx + sin V3 + 7 = 0, (7a) 
V5z(0) = (Px{L) + fisiinp{L) = he- (7b) 

We note that independently of the junction geometry the presence of a point inhomoge- 
neity at the right end leads to nonlinear boundary conditions. 

Another example of an inhomogeneity that will be considered below is a narrow rect- 
angular well of width A <C L and depth 1 — k: 

3d{x) = l-{l-t,)e{x + A-L). 

Here 9{x) is the Heaviside step function. The parameter k > specifies the fraction of 
the Josephson current through the inhomogeneity. We note that the inhomogeneity of 
the barrier layer represents a microresistor only for k < 1. In particular, if k = the 
Josephson current is absent on the segment [L — A, L]. The choice k = 1 means that the 
thickness of the barrier layer of the junction is constant, and the amplitude joix) = 1 for 
X G (0, L]. For k > 1 the inhomogeneity is a microshort that repels magnetic flux vortices 

The corresponding boundary-value problem for the magnetic flux ip{x) in the case of 
overlapping geometry is written in the form 

->fxx+jD{x)sm(f-'j = 0, (8a) 
<^x{0) = he, iPx{L) = he. (8b) 

In a numerical simulation of the magnetic-fleld dependence of the critical current, 
transitions from the superconducting to the resistive regime are interpreted mathemati- 
cally as bifurcations of the magnetic flux in the junction upon variation of the parameters 
P^ . For studying the stability, each solution (p{x) of the boundary- value problems (jH)) 
is associated to a regular Sturm-Liouville (S-L) problem with a potential that depends on 
the concrete solution (see the details in Refs. [11] - [IB])- The S-L eigenvalue problems 
allow one to assess the stability or instability of (p{x). The boundary conditions are 
determined by the choice of model. In particular, for model ((Tj) with a point inhomogeneity 
the S-L problem takes the form 

-'ipxx + q{x) ij = X^J , (9a) 
^,(0) = 0, V^,(L)+/icos¥^(L)7A(L) = 0, (9b) 
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(10a) 
(10b) 



^\x) dx-l = 0. 



(11) 







We note that problems and (fTUjl have a discrete non-degenerate spectrum bounded 
from below. In addition, since the solutions of the boundary-value problems ((Zj) and © 
depend smoothly on the parameters Hb and 7, the eigenvalues and eigenfunctions of 
equations and (fTUj) also depend on these parameters. Let Amin(/iB,7) be the minimum 
eigenvalue corresponding to a certain distribution v^(x). The point with coordinates hs 
and 7 on the {hB,'y) plane is called a bifurcation point of ipi^x) if the following 
condition holds: 



The bifurcation curve of any magnetic flux distribution in a junction is the geometric 
locus of bifurcation points of a concrete vortex, corresponding to pairs {Jib-,'^) for which 
Eq. dig) holds. 

In view of the fact that for specified Kb and 7 the boundary- value problems (|7j) and 
dnj, as a rule, have more than one solution, the critical curve of the junction as a whole is 
constructed as the envelope of bifurcation curves for different solutions. In other words, 
the critical curve of the junction consists of pieces of the bifurcation curves corresponding 
to different distributions with the highest critical current 7 for the given field Kb- 

For calculation of relation (fT^ we used the method proposed in Ref. [T7]: the system 
of equations ((Tj) and © [respectively Eq. (jH}, (fTUj) ] at fixed A and field Kb (and current 
7) is solved as a nonlinear eigenvalue problem with a spectral parameter 7 (respectively, 
Kb)- This method of seeking bifurcation points and its generalization was used in Ref. 
[TH] for solving a wide class of physical problems. 



Let us compare the results of numerical simulations for constructing the critical curves 
for the model of a JJ in an overlapping geometry with an exponentially varying width 
and a model with an inhomogeneity on the right end. 

FigureHshows the "bifurcation" distributions of the magnetic field ^^^^(x) for the single 
fiuxon $^ along the junction for three JJ models (EJJ, IJJ, and 5JJ) for Amin = 10~^, i.e., 
before the destruction of the fiuxon by the current 7. The corresponding distributions of 
the Josephson current joix) smip in the J J are presented in Fig. |21 

We note the good qualitative agreement of the curves. Their quantitative difference 
is observed predominantly in the region where the amplitude of the Josephson current is 
minimum; this is due to the form of the inhomogeneity (attractor). 



Amm(7, ^e) = . 



(12) 



4. NUMERICAL RESULTS 
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Figure 1: The magnetic field of tlie fluxon 
$^ for three models of the inhomogeneity. 



Figure 2: Distribution of the Josephson 
current for three models of the inhomo- 
geneity. 



Figure El demonstrates the critical curves of the form (fT^ corresponding to stable 
solutions of the boundary- value problems and (jH))- The solid curves correspond to 
distributions of the magnetic flux for the model of a JJ with a width that varies by 
an exponential law (EJJ) with a coefficient a = 0.07, the dashed curves are for models 
of an inhomogeneous JJ with a finite rectangular inhomogeneity (UJ) with parameters 
A = 0.7, K = 0.1. Each critical curve refiects a certain closed figure on the (/tB,7) plane, 
the currents hs = /imin and hs = /imax are the lower and upper critical magnetic fields for 
the vortex under study. We note that because of the presence of an additional "geometric" 
current cr(v9^ — hs) in the EJJ model the critical value of the current 7 corresponding to 
the two indicated values of hs is nonzero. 
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Overlap geometry: L = 7 

EJJ: a = 0.07 
IJJ: A = 0.7, K = 0.1 
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Figure 3: Bifurcation curves for a JJ of variable width and a JJ with a rectangular 
inhomogeneity. 



A numerical simulation shows that, unlike a J J with a resistive inhomogeneity within 
the segment (0, L) (Refs. [T7j and [H]), in a junction with an inhomogeneity at the end, 
even for rather large values of (the solutions were checked numerically to hs = 10) 
there are no stable mixed fiuxon-antifiuxon vortices contributing to the critical curve. 
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Therefore the critical curves of the junction as a whole consist only of segments of the 
critical curves for a Meissner (M) distribution and "pure" n-fluxon vortices n = 
1, 2, . . .. The maxima of the critical curves fall off monotonically with increasing field hs 
(in junctions with an internal inhomogeneity the maxima of the critical curves fall off 
non-monotonically because of the stabilizing influence ^H] of the inhomogeneity on the 
mixed fiuxon-antifiuxon vortex pairs). 

Let us illustrate the process of vortex destruction (the transition to an unstable state) 
upon variation of the current 7 for the case of the IJJ model. In a field /i^ = 1 at a 
current 7 > 0, two stable distributions of magnetic flux can exist in the model: M and 
As the current increases, the main fiuxon $^ is destroyed first (the critical current for 
it is 7cr($^) ~ 0.062), and then the Meissner distribution. Thus for Hb = 1 the critical 
current of the junction 7^ = 7cr(M) ~ 0.597. For hs = 1.9 the boundary-value problem 
(jS} has three stable solutions — M, and which occur in the following succession 
with increasing 7: $^ ^ A<f — > The critical current of the junction in this case is 
7cr = 7cr('^'^) ~ 0.157. Analogously, for Hb = 2.2 we have 7cr7cr('^'^) ~ 0.092. 



Figure 4: Critical curves of the current versus magnetic field for three model inhomo- 
geneities. 

The critical curves of the current versus magnetic field for the three JJ models discussed 
above are exhibited in Fig. |3J The dot-and-dash line corresponds to the model of a JJ 
with a point inhomogeneity (/i = 0.2) on the right end (5JJ). Note the good qualitative 
agreement of the critical curves for the three models. For values of Hb that are not too 
high, the leftward displacement of the critical curves of the individual vortices by the 
"geometric" current in the EJJ model has only a weak influence on the behavior of the 
critical curve of the junction whole. 

We note that variation of the parameters of a rectangular inhomogeneity (of width 
A and depth n) has a weak influence on the numerical results, i.e., the boundary- value 
problem (jSJ is structurally stable upon variation of those two parameters. This assertion 
is illustrated in Fig. |31 which shows the dependence (fT^ for the main fiuxon $^ at 7cr > 0, 
K = 0.1 and values of A = 0.5 and A = 1. Analogously, Fig. El shows the dependence (fT^ 
for the main fiuxon $^ for ^ 0, A = 0.7 and parameter values k, = 0.1 and k, = 0.8. 
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Figure 5: Influence of the parameter A on Figure 6: Influence of the parameter k on 



In this paper we have shown that every JJ of variable width can be placed in correspon- 
dence with a quasi-one- dimensional junction with a variable thickness of the barrier layer. 
In the case of junctions whose width varies by an exponential law, the barrier layer of 
the corresponding quasi-one-dimensional JJ contains a resistive inhomogeneity distributed 
along the junction and acting as an attractor for magnetic flux vortices. 

The curves of the critical current versus magnetic field obtained from the numerical 
simulation demonstrate that such an inhomogeneity can be replaced by a local inhomoge- 
neity at the end of the JJ, which may have certain advantages from a technological point 
of view. 

The authors thank Prof. Yu. Kolesnichenko (B. Verkin Institute for Low Temperature 
Physics and Engineering, Kharkov) for helpful discussions. 
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